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We show that (~ l/(Q 2 ) fc ) power corrections to the spin structure function gi at small x and large 
£ — . Q 2 are generated perturbatively from the regulated infrared divergencies. When Q 2 are small, the 

corrections are ~ (Q 2 ) k ■ We present the explicit series of such terms as well as the formulae for their 
resummation. These contributions are not included in the standard analysis of the experimental 
. data. We argue that accounting for such terms can sizably change the impact of the power corrections 

( — I ■ conventionally attributed to the higher twists. 



3 

00 



o 

O 



43 



X 



PACS numbers: 12.38.Cy 



I. INTRODUCTION 



(N 
> 

The theoretical description of the Q 2 -dependence of the structure function g\ in perturbative QCD is mostly 
performed for the kinematic region of large Q 2 . However, for the phenomenological analysis of the results of the 
COMPASS collaboration one needs the possible knowledge of g% at small x and small Q 2 (see e.g. Ref. [l[). Let 
us remind that the conventional Standard Approach (SA) based on combining the DGLAP evolution equations^ 
with phenomenological inputs for the initial parton densities [H, cannot be used for description of g\ in this region. 
Strictly speaking, the SA can be applied only for large values of x (x ~ 1) and large Q 2 : Q 2 ^> Q 2 , where Q\ is the 
starting point of the Q 2 -evolution. Indeed the small- a: region lies beyond the reach of the SA because DGLAP does 
not include the total resummation of oc ln fc 1/x terms. In order to describe the experimental data at small x, in the 
SA one has to include singular oc x~ a terms in the expressions for the initial parton densities. Such factors act as 
the leading singularities (simple poles) in the Mellin transform of gx(x) and provide g\ with the Regge asymptotics: 
gi ~ x~ a when x — ► 0. 

On the other hand, in our approach [J Q based on the total resummation of the leading ln fe 1/x terms in perturbative 
QCD, the Regge behaviour of g\ at x — ► appears naturally, independent of the value of Q 2 . Indeed it results from 
the leading singularities of the anomalous dimensions and coefficient functions, which are branching points and not 
simple poles. The singularities of the anomalous dimensions and coefficient functions coincide. This is very important 
because it guarantees the independence of Q 2 of the intercepts of g\. Furthermore by fitting the experimental data 
in our approach one does not need ad hoc singular factors ~ x~ a in the initial parton densities. 

The theoretical aspects of the power corrections to the DIS structure functions were considered in Refs. @ 
at large x and Q 2 . The interplay between the perturbative and non-perturbative corrections in the hard kinematics 
was recently considered in detail in Ref [9(. In fitting experimental data in the small- 2 region, the discrepancies from 
the SA predictions are conventionally interpreted (see e.g. Refs. [l(|) as the (non-perturbative) higher twist power 
(~ (1/Q 2 ) k ) contributions. However as the small- x region is beyond the reach of SA, the real size of the higher twist 
corrections can be erroneously overestimated. It is clear that systematic account of power corrections ~ (1/Q 2 ) fc can 
be satisfactory only when using formulas that already account for the resummation of logarithmic contributions. 

We have recently proposed in Ref. [ll| a generalization of our previous results 0,01 for gx. Although this extension 
goes beyond the logarithmic accuracy we keep, it looks quite natural. In addition to the standard region of large Q 2 , 
this generalization can describe the small-Q 2 region, though in a model-dependent way. Our suggestion, based on the 
analysis of the Feynman graphs for g\ at small x, is to replace Q 2 by Q 2 + /i 2 in our previous formulas, with /i being 
the infrared cut-off Such a shift also leads to replacement the variable x byx' = x + /j 2 /2pq. The variable x' is similar 
to the Nachtmann variable 1 . 



We are grateful to L.N. Lipatov for reminding this. 
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In the present paper we show that regulating the infrared divergencies originates also the power corrections for gi 
at small x. They differ from the well-known power corrections [E 0, H| related to the resummation of the Sudakov 
logarithms: first, they come from the ladder Feynman graphs in the Regge kinematics where virtual gluons are not 
always soft; then, in contrast to Refs. [1, 0> fl @ j we never use 2 the parametrization a s = a(Q 2 ) and do not consider 
inclusions of power terms (l3j into the standard expressions for a s . We present the explicit formulas for the total 
resummation of these new corrections and show that at large Q 2 they are ~ l/(<5 2 ) fc but in the small- Q 2 region they 
are ~ (Q 2 ) k - The paper is organized as follows: in Sect. 2 we explain why there should be the infrared perturbative 
and non-perturbative power corrections to g\. We also discuss in this Sect, the difference in the infrared properties 
of gi at large and small x. In Sect. 3 we remind the essence of the model for gi suggested in Ref. [ll[ and list explicit 
expressions for the singlet and non-singlet gi at the kinematic region of small x and arbitrary Q 2 . In Sect. 4 we give 
in more detail than in Ref. fTlj theoretical arguments in favor of those expressions. The expressions for gi enlisted in 
Sect. 3 and proved in Sect. 4 account for the total resummation of leading logarithms of x and Q 2 when Q 2 are large 
but they are also valid for small Q 2 . They implicitly include the power Q 2 -corrections. The explicit expressions for 
the power corrections are extracted from those formulas in Sect. 5. We show here that the power corrections at large 
and small Q 2 are quite different but the lowest twist contribution to gi is always leading regardless of Q 2 . Finally, 
Sect. 6 is for discussion. 



II. ORIGIN OF THE INFRARED-DEPENDENT CONTRIBUTIONS TO gi 

It is well-known that the DIS structure functions are introduced through the hadronic tensor Wy. v . In particular, 
the spin-dependent part of W^ v for the electron-proton DIS is parameterized by the structure functions g\ and g2- 

S p9l (x,Q 2 ) + (s p ~^) g 2 (x,Q 2 )] (1) 
V pq / J 

where p, M, S are the proton momentum, mass and spin respectively, q is the virtual photon momentum. The spin 
structure functions gi.2 have the non-singlet, g±§ ■> and singlet, gf 2 , components. In general, gi^ (and all other DIS 
structure functions) acquire both perturbative and non-perturbative QCD contributions. 

In the first place, there is a non-perturbative term W^ U PT in the region of small Q 2 . For instance, there are known 
examples of the lattice calculations for some structure functions (see e.g. Ref. [Hi), although basically W pl , PT is 
poorly known. 

The other part, W p J ' , of W^u includes both perturbative and non-perturbative QCD contributions. The standard 
way to calculate the structure functions is using the factorization. According to it, W pv is regarded as a convolution 

W?J = W*, ®% + W° v ® $ 3 (2) 

of the partonic tensors WZ V , W® v (where q, g label the incoming parton, i.e. a quark or a gluon respectively) and 
probabilities $ 9 £( to find this parton (quark or gluon) in the incoming hadron. The probabilities $ q , g include both 
perturbative and non-perturbative contributions whereas purely perturbative tensors W^]f describe their x- and Q 2 - 
evolutions. There are no known explicit expressions for $ g , g . Instead, they are approximated by the initial parton 
densities 5q and Sg defined from fitting experimental data at Q 2 ~ 1 GeV 2 and x ~ 1. The partonic tensors 
W^ v , WR V evolve these densities into the region where Q 2 ^> fi 2 and x <C 1, with /i 2 being the starting point of the 
Q 2 -evolution. Such a contribution to W pv is called the lowest twist contribution W p J '. Besides, there are the higher 
twists contributions W^J . They can be interpreted either in terms of more involved convolution or as essentially 
non-perturbative objects. 

In order to calculate W^, one has to regulate infrared (IR) singularities in involved Feynman graphs. Such a 
regulation is different for large and small x. At large x one can use DGLAP. In DGLAP this problem is solved with 
assuming a non-zero virtuality /i 2 for the initial partons and imposing the ordering 

M 2 < k\ ± < k 2 <k 2 n± <Q 2 (3) 

on the transverse momenta of the ladder partons (the numeration in Eq. ([3]) runs from the bottom to the top of the 
ladders). Eq. §3§ manifests that k r j_ acts as an infrared cut-off for integrating over k r+ \. However, the ordering allows 
to collect the contributions that are leading at large x only. In order to account for the leading (double-logarithmic) 



2 It is known (see Ref. |12|| ) that this parameterizations fails for gi at small x 
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contributions at x <C 1, the upper limit of the integration in Eq. ([3]) should be changed for (p + q) 2 « 2pq and the 
ordering should be lifted. Without the ordering, the IR singularities in the ladder rungs are not any longer regulated 
automatically, so an infrared cut-off should be introduced for integration over every loop momentum, ft is the reason 
why the /z -dependence of g\ is getting more involved at small x. Usually the cut-off is identified with /z, the starting 
point of the Q 2 -evolution, though it is not mandatory. Obviously, the value of /z should be large enough to justify 
using the Perturbative QCD: 



*i >l?> &QCD- (4) 

Generally speaking, there are different ways to introduce infrared cut-offs but providing IR-divergent propagators 
with fictiones masses is most wide-spread. Obviously, no observables should depend on a value of /z and ways of its 
introducing. It means that the explicit /z -dependence in W5j? should be compensated by a /z -dependence of Sq, Sg. 
However, the latter are known as phenomenological expressions containing a set of numerical parameters fixed from 
fitting experimental data. Those parameters are supposed to be fi -dependent, though in unexplicit way. 

Before considering the IR properties of g\ at small x, let us discuss the well-known expression for the non-singlet 
<7i in the Standard Approach: 

r°° rlii) r f® 2 dk 2 i 

g? s (x,Q 2 7 n 2 ) = (e 2 q /2)j^ 2L x -»C t (u)5q(L,) exp [j ^ ^7 m (k 2 x , a s (kl))\ (5) 

where e q is the quark electric charge, C q is the coefficient function and "f qq is the DGLAP non-singlet anomalous 
dimension. When ^ qq is taken in LO, the integral in Eq. ([5]) is known to be (1/6) m[ln(Q 2 //i 2 )/ ln(fi 2 /Aq CD )], with 
6 being the first coefficient of the Gell-Mann -Low function. The expression for the singlet g\ looks similar, though 
more involved. Both the coefficient functions and the anomalous dimensions in SA are known in few first orders 
in a s . It means that using Eq. ^ is theoretically based for the kinematic region where x are not far from 1 and 
Q 2 3> /z 2 ~ 1 GeV 2 . The expression for g% in Eq. ([SJ depends on the value of /z. This dependence is supposed to 
disappear when g\ is complemented by a contribution g^ T extracted from W^J '. In practice, the treating experimental 
data on Polarized DIS is carried out as follows (see e.g. the recent paper Ref. [lCj] and refs therein): the data are 
compared with g\ of Eq. |(5J) and the discrepancy is attributed to the higher twists impact. However, SA for g\ is 
reliable at large x. At the small- a; region it should be modified. In the next section we present explicit expressions 
for gf and g± s replacing the DGLAP expressions in the small- x region. 



III. EXPRESSIONS FOR gi AT SMALL x AND ARBITRARY Q 2 

When x -C 1, the contributions ~ hi (1/x) are large, so they should be accounted for to all orders in the QCD 
coupling. The total resummation of the leading logarithms of x for gi was done in Refs. for the region of Q 2 > /z 2 . 
We remind that, contrary to DGLAP, the expressions for g\ in Refs. [1, [f| are valid both for large Q 2 and for Q 2 ~ /i 2 . 
Recently, in Ref. [ll| we have suggested a simple prescription to generalize those results to arbitrary values of Q 2 : 
in the formulas of Refs. [1, H[ Q 2 should be replaced by Q 2 + is 2 . This conclusion follows from the observation that 
the contributions of Feynman graphs to g\ at small x depend on Q 2 through Q 2 + fi 2 only. It automatically leads to 
the shift x — > x + z, with z = /z 2 /2pq. As the prescription is beyond the logarithmic accuracy that we kept in our 
previous papers, we call it a model. The theoretical grounds of this model arc given in the next Sect. According to 
our results, the non-singlet, g^ s ', component of g\ at the small- x region is given by the following expression: 

+ + " 2 » - <^> £ ^y c ^<^r M ■ <*> 

with the new coefficient functions C^vs, 



Cns(u) = ^—^ (7) 

LU - H ns {uj) 



and the anomalous dimensions Hns, 



H NS = (1/2) u - V^ 2 ~ B{uj) (8) 
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where 



B{uj) = (4ttC f (1 +lu/2)A(lu) + D{uj))/(2tt 2 ) 



(9) 



D(u) and A(ui) in Eq. ([9j> are expressed in terms of p = ln(l/x), r\ = \n(p 2 / Aq CD ), b = (33 — 2n/)/127r and the color 
factors C F = 4/3, AT = 3: 



D{u>) = 



2C F 



A(u) 



P + V 



1 



^2 _|_ 



(p + v) 2 + 

dpe^P 



P + V 



[p + vf 



(10) 



(11) 



-Hjvs and Cjvs account for DL and SL contributions to all orders in a s and, contrary to the DGLAP phenomenology, 
Sq does not contain singular factors. The infrared cut-off p obeys Eq. Expression © is valid for large Q 2 , i.e. 
for Q 2 ^> p 2 where x ^> z and for small Q 2 , Q 2 < p 2 where x < z. The expression for the singlet component, gf , of 
gi is more involved: 



g? = g[ +) + g[-\ 



(12) 



with 



<e q > 



1 



2m \z + x 



A' 



q q 2nuo 2 9 



Q 2 + p 
p 2 



(13) 



where < e 2 > stands for the sum of electric charges: < e 2 >= 10/9 for n/ = 4, Sq and Sg are the initial quark and 
gluon densities. 

The exponents fi(±) and coefficient functions Cqjg are: 



n (±) = ^ [H m + H gg ± R] 



c« = 


UJ 

RT 


(Hgg 


-«(-))( 




UJ 

rT 




-Hgg){ 


C (+) = 


UJ 

RT 




-«(-))( 


cV = 


UJ 

RT 







H gg)i u - H qq) - H qg H gq + H qg(^(+) - w ) 



.4' 



2ttw 2 

A' 
2vrw 2 



(14) 
(15) 



Here 



7? 



[H qq -H gg f + AH qg H g 



T - UJ 2 ~ Uj(Hgg + Hgg) + (HggHgq - HgqHqgJ 



(16) 



and 



H gg = - 



w + z + bgg 6gg 



w + z - bqq b " 



tt _ u qg 

qg ~ ~Y 

TT - bgq 



(17) 



where 



Z=-J=\I(UJ 2 - 2{b qq + b gg )) + y/(<J» - 2(bgg+b gg )) 2 - A(bgg - b g g ) 2 - Wbggbgg 



(18) 
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with 



bik — O-ik + Vik 



The Born contributions a,i- are defined as follows: 



A(u)C f 



a l<i ~ r. 
At last, non-ladder contributions are: 



2ir 



A'{uj)C f 



"gg 



<<HJ 



ANA(co) 
2^ 



with 



'99 



Cf_ 
2N 



m gg = -2N Z 



D(w) 



N 



'99 



'</y 



= -NC f 



The additional factor ^ — j in the coefficients in Eqs. (fT5)) . with 



A'M = - 



P 



o G° + ^) 5 



(19) 
(20) 
(21) 
(22) 

(23) 



is the small-w estimate for the quark box diagram that dominates in the Born term relating initial gluons to the 
electromagnetic current. A'(u>) is the Mellin representation of the QCD running coupling a s involved in the quark 
box. 

Besides resummation of the leading logarithms, Eqs. (|6ll2p differ from DGLAP in the parametrization of a s : the 
DGLAP -prescription is a s = a s (Q 2 ) whereas in our approach a s is replaced by A and A' defined in Eqs. (|11I23[) . 
Such a difference results into a drastic difference in the form of the Q 2 -dependence of g\ between our approach and 
SA: instead of the factor (Q 2 + fi 2 )/fi 2 in Eqs. (jSTT2"l) , the SA leads to ln(Q 2 / A^ c D ) . 



IV. THEORETICAL GROUNDS FOR THE SHIFT Q 2 -> Q 2 + fi 2 IN EQS. © AND (f!2\i 

Both the singlet and non-singlet components of gi obey the following Bethe-Salpeter equation depicted in Fig. 1: 

.91 - s? orn + * J j^yA~^) 5 ^ + fc ) 2 ~ m ^ {k 2_ ml? M ^ fc ) ( 24 ) 

where the 8 -function (together with the factor —2m) corresponds to the cut propagator of the uppermost quark with 
momentum k and mass m q coupled to the virtual photon lines and 2k\ appear after simplifying the spin structure 
of the equation. M(p, k) stands for the cut parton ladders and the initial parton densities. In other words, M(p, k) 
incorporates both the initial parton densities and radiative corrections to them. This object can be called the polarized 
parton distribution function. In the present paper we will address it simply as the (parton) distribution, skipping 
other words. For the sake of simplicity we dropped unessential numerical factors (e 2 /2 for f NS and < e 2 > /2 for the 
singlet) in Eq. Obviously, M in Eq. (pM)) cannot depend on Q 2 . 



A. Prescription for the IR regularization of M 



The IR-divergent contributions to M should be regulated. We follow the standard prescription and assign a fictitious 
mass [i to the gluons in the IR-divergent propagators. In particular, in the ladder graphs such propagators are the 
vertical ones. We assume that the value of fi satisfies Eq. (|4|). In contrast to the gluon ladders, quark ladders are 
IR-stable because the quark mass m q acts as an IR cut-off. In order to use the same cut-off /x for regulating both 
gluon and quark IR divergences, we assume that, in addition to Eq. ((4]), fi 3> m q . After that m q can be dropped. 
Therefore in order to regulate IR singularities, we should insert /i 2 in the IR -divergent, strut (vertical) propagators, 
for both quarks and gluons. The horizontal propagators (rungs) are IR-stable. This converts Eq. (|2~1| into 

/r1 4 k 9k 2 
-^(-2m)6((q + fc) 2 ) _^ )2 M(pk, k 2 + fi 2 ) (25) 
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B. Solving the Bethe-Salpeter equation (|25[) 

As the kinematics x *C 1 is of the Regge type, we need to express M in the Regge kinematics as well. Let us notice 3 
that the expressions for M can be obtained from our results for gi in Ref. 0, Q with replacing the external photon 
virtuality q 2 = —Q 2 by the external quark virtuality k 2 and x by —k 2 /wa. In the first place, we focus on applying 
Eq. {25 to gf s , the non-singlet part of g\ and denote M NS the involved quark distribution. The expression for M NS 
accounting for the total rcsummation of leading logarithmic contributions can also be borrowed from our formula for 
Qi s obtained in Ref. [f|. This expression reads: 

duJ ( 2pk \w x /-fc 2 + /i 2 \Hns(«j) 



/DC 



too 



M NS (p,k)= I ^[^=f— 2 ) uf( U )6q{ U )l— -^-) (26) 



where / = 8tt 2 Hns, 5q{u>) is the initial quark density in the u -space and Hns is given by Eq. ([8]). H^s an d / include 
the total resummation of leading logarithmic contributions. Similarly, the expression for the singlet distribution M 
can be obtained from our results for gf in Ref. [4J. It is convenient to rewrite Eq. (p?5|) in terms of the Sudakov 
variables for momentum k: 

k = -aq' + (3p + k ± , (27) 

with q' = q + xp, so that q' 2 « p 2 » 0. Substituting Eq. (f2"6"| into Eq. (|23|) . using the Sudakov variables and changing 
the order of the integrations, we obtain the following equation for gi S : 

-| /"ZOO 7 f 7,2 

= + _ j _ uf ( u ) Sq ( u) J dadf3 dk 2 ± (wQ/? + fc L 2 +M 2 ) 2 ^M + ^« " W( *P ~ k l- C 28 ) 

/ wa \ w /waft + k\ + /j, 2 \ h ns 
\wa/3 + kj_+ n 2 ) \ fi 2 ) 

where we have denoted w — 2pq. As we consider i < 1, we can neglect xa compared to (3. Using the 5 -function for 
integration over /?, we arrive at 

NS Born 1 P°° duj f dadk\ / wa \" ( aQ 2 + k 2 , + fi 2 \ h ns 

The integration region in Eq. (|25p is shown in Fig. 2. It is outlined by the following restrictions: 

w > fc 2 > 0, aQ 2 + k]_ + fi 2 . (30) 

Integrating over a and fc 2 ^ yields different contributions, depending on the ratio between aQ 2 and k\. The most 
important contribution comes from the region D in Fig. 2. After integration over a in this region we get 

9i =9i — w/ w <5 ? w - / 2 - x f 2 1 ( 5 — (31) 



.Born 



l± 2 \kj_ + [I 2 ) V £i 2 
1 r°° dbJ ,„ , , /" U ' +m2 dtfW\»f t \Hns 



9r ' " + ^ i_ 2^ mSq{uj) t(t) (?) ^ 

The leading contribution in Eq. (fBTj) comes from the lowest limit t — Q 2 + fj, 2 and gives 

,32) 



» s _ Bm , i r°° <<" /(»)*) » rr<5 2 + f^"" 

yi — ai 



8^y_ loo 2^(a.-^5)V(g 2 + A i 2 )^ v M 2 

which proves the validity of the shift Q 2 — > Q 2 + /i 2 suggested in Ref. [ll|. Let us brush up Eq. ([32l . Replacing /(w) 
by 8tt 2 H ns we see that in Eq. (|3"2")l 

-^-=-1 + — (33) 

The first term in Eq. (|33|) cancels the Born contribution g® orn and the second term is the non-singlet coefficient 
function (see Eq. ©). Therefore, we arrive at Eq. ^ for the non-singlet g^ s at small x and arbitrary Q 2 . Eq. (TT 
for the singlet g\ in the same kinematic can be proved similarly. 



3 We will consider such distributions in more detail in next our paper. 
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V. INFRARED POWER CORRECTIONS AT SMALL x 
A. Power corrections at large Q 2 

In the kinematics where Q 2 > /i 2 and therefore x > z, the terms with Q 2 + fi 2 in Eqs. (|6I12[) can be expanded into 
series in fj?/Q 2 : 



1 



x + z 



) 2 + [1 2 \ h ns 



\x + z) \ /i 2 / \x) V /j, 2 J 



f \ Hj 

-Q 2 \"± 



(i)"(l)""hi;^K| 



(34) 



fe=l 



where 



— 



(-L0 + H NS )(-U + ffjys - + gjV5 - k + 1) 

fc! 



(35) 



T, 



( ±) _ (-a; + n ± )(-uj + n±- + n ± -k + i) 

fc! 



It allows to rewrite Eqs. (|6I12|) as follows: 

g? s (x + z,Cf) = ~g? s {x,Q 2 ) + ~g^{x/y,Q 2 )®Y,^ 2 IQ 2 ) k E k S (v), 



(36) 



fc=i 



(x + z 7 Q 2 ) = gf(x,Q 2 ) + ^( M 2 /Q 2 ) fc [?^ +) («/!/, Q'JsE^^+^fir/tf.Q 2 )®^^) 



k=l 
-,NS OT1 j ~ - 



where, using the conventional terms, g{ and gf can be named the non-singlet and singlet components of the lowest 
twist contribution to g\\ 



g? s 



9i 



(e 2 q /2) 



du) (1 
2m \x 



(-Yc N s(Lj)Sq(u;)(- 



Q2 > ffjvs ( u ) 



(37) 



g[ +) - 

A' 

27TW 2 



«1 — 



2 

' Q 2 N ° 
/i 



du) /I 
2m \ x 
2 N n 



) 2 \ n 



and 



(38) 



The right-hand sides in Eq. (|36|) are the products of the power corrections and g±. We call these corrections the 
infrared power corrections. The functions gf , g^ s in Eq. (|37[) were obtained in Refs. [1, H| and they correspond to 
the lowest twist contribution. They differ from the lowest twist DGLAP expressions for gi by the total resummation 
of the leading logarithms of x and by the new parametrization of a s given by Eqss. (| 1 1 123[) . They include the most 
important at small x contributions of the LO and NLO DGLAP formulas. When the lowest twist expressions of 
Eq. (|37[) are used for analysis of experimental data of the polarized DIS, the power series in the rhs of Eq. (|3"6"|) look as 
new independent contributions. However, the left-hand sides of Eq. (13"6")) account for the total resummation of these 
corrections. Finally, let us notice that the infrared power corrections in Eq. (|36[) have nothing to do with the standard 
parametrization a s — a s (Q 2 ) as we do not use it. 
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B. Power corrections at small Q 2 



When Q 2 < /i 2 , g^ s and gf cannot be expanded similarly to Eq. (|36[) . The power corrections for small Q 2 are 
different. Indeed, in this case 



\x + z 



Q 2 + H 2 \ H *s ( \ 

z 



ix + z) ( 



fc=i 



fc=i 



2 X fe 



(39) 



ft leads to the following expressions for gi at small Q 2 



g? s (x + z,Q 2 ) = ^(z^VsfWw^EtW^ffe)- 



(40) 



fc=i 

ffl s (x + z,Q 2 ) = ^W^ + ^Q 2 //* 2 )* [^(^^(gl^^+^r^Vtf,/* 2 )®^"^!/) 



where the lowest twist contributions g^ s , gf are given by Eq. (|3"T| and , are defined in Eq. ([55)1. Both 
and gf do not depend on a; and Q 2 . Instead, they depend on the total energy (p + q) 2 of the process and are constants 
when the 2pq is fixed. Both the x and Q 2 -dependence are now associated with the power corrections. Eqs. (|36l40p 
show that the infrared Q 2 - corrections are different for large and small Q 2 , so that the series of Eq. (l40l cannot be 



extrapolated into the region of small Q 2 and similarly the series in Eq. cannot be extrapolated into the large Q 2 - 
region. Besides the terms with {Q 2 ) k given by Eq. (|40[) . similar contributions can come from other sources which are 
beyond our control. However, the {Q 2 ) k -terms in Eq. (|4U|) are multiplied by the functions g^ s and gf which include 
the total resummation of the leading logarithms and therefore they are supposed to dominate, at small x, over the 
coefficients at the other (Q 2 ) k -terms. 

Although the large-Q 2 expansion (f3"4"| and the small-Q 2 expansion (13"9"P look quite similar, the power corrections to 
gi s are actually different for large and small Q 2 . It is easy to check that the linear in /j, 2 /Q 2 term is present in the 
large-Q 2 expansion of Eq. (|36|) for g^ s while the term with Q 2 //j, 2 is absent in the gf^-expansion of Eq. (|40| . 



VI. DISCUSSION 



At the region of small x, the DGLAP ordering §5§ should be lifted for accounting for leading logarithms of 1/x, 
so the infrared cut-off /j, should be introduced explicitly to regulate the infrared divergencies in every rung of the 
Feynman graphs contributing to g\. Similarly to DGLAP, /i can also play the role of the starting point of the Q 2 
-evolution, though not obligatory. With both the perturbative and non-perturbative contributions accounted for, g\ 
does not depend on fi. However, this dependence does exist when the non-perturbative contributions are neglected 
or accounted for implicitly through the fits for the initial parton densities. In this case the structure function g\ 
depends on the value of \i and the way it has been introduced. Introducing \i as the fictitious mass inserted into the 
IR -divergent propagators, leads to Eqs. J6II2[) suggested in Ref. [ll| for g\. The expressions of Eqs. (|6lf 2\ include 
the total resummation of double-logarithms and the most important part of single-logarithms of x. They are obtained 
from our previous results with the shift Q 2 — > Q 2 +/i 2 . The theoretical grounds for such a shift are given by Eqs. (|9l32p. 
Eqs. (|6|12p . in contrast to DGLAP, can be used both for large and small Q 2 . Having been expanded into the series in 
l/(Q 2 ) k at large Q 2 (or into series in Q 2 at small Q 2 ), Eqs. (|6ll2p yield the power Q 2 - corrections. The series of the 
corrections are represented by expressions (I36I40[) . The power series of Eqs. Q36I40P for large and small Q 2 are derived 
from the same formulas. However after the expansion has been made, they cannot be related to each other with 
simply varying Q 2 . We suggest that accounting for the new source of the power contributions given by Eqs. (|6ll2p 
can sizably change the conventional analysis of the higher twists contributions to the Polarized DIS because such 
contributions appear in the present analysis of experimental data as a discrepancy between the experimental data and 
the Standard Approach predictions. Let us remind that in Ref. [l5[ we showed that the singular (~ x~ a ) factors in the 
standard fits mimic the total resummations of ln fe x, i.e. they have a purely perturbative origin contrary to previous 
common expectations. Similarly, a good portion of the commonly believed non-perturbative power corrections in 
the conventional analysis of experimental data can actually be of the perturbative infrared origin. However being 
misinterpreted as non-perturbative terms, they can mimic the power expansion in Eq. (|36p . In particular, Eq. (JBJ) 
predicts that the power ~ l(Q 2 ) k -corrections to g^ s should appear at Q 2 > 1 GeV 2 and cannot appear at smaller 
values of Q 2 . It agrees with the phenomenological observations obtained in Refs. [To| from conventional analysis of 
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experimental data. On the other hand, Eq. ([12]) predicts that the similar power corrections to the singlet g\ should 
be seen at greater values of Q 2 . Clearly, the use of Eqs. (|6I12[) for the lower twist contributions to <? 1; instead of 
DGLAP, would allow one to revise the impact of the genuine higher twists contributions which are known to be of 
the non-perturbative origin. 

Finally, we would like to remind that our results explicitly depend on the infrared cut-off /i. As pointed out in 
Sect. II, such a dependence would vanish if analytic expressions for the probabilities $> q ,g were obtained and used in 
Eqs. (|6ll2p instead of Sq and 8g. 
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2: The integration region over a and k\ in Eq. 



